Unified Treatment of Even and Odd Anharmonic Oscillators of Arbitrary Degree 
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We present a unified treatment, including higher-order corrections, of anharmonic oscillators of arbitrary even 
and odd degree. Our approach is based on a dispersion relation which takes advantage of the "PT-symmetry 
of odd potentials for imaginary coupling parameter, and of generalized quantization conditions which take into 
account instanton contributions. We find a number of explicit new results, including the general behaviour of 
large-order perturbation theory for arbitrary levels of odd anharmonic oscillators, and subleading corrections to 
the decay width of excited states for odd potentials, which are numerically significant. 
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Introduction. — One-dimensional anharmonic oscillators 
are quite basic. Because of their enormous phenomenologi- 
cal significance, they occupy a unique position within quan- 
tum theory JT} [2j [3j |4l [5J . They are treated at various levels of 
sophistication in nearly every textbook on quantum mechan- 
ics. Here, we use the Hamiltonians of even oscillators in the 
convention 

H N (g) = -\d 2 q + \q 2 + gq N , (TV even) , (1) 

and odd oscillators as follows, 

h M {g) = -\d 2 q + \q 2 + ^q M . (M odd) . (2) 

For g < 0, the potential of an even oscillator has a double- 
hump structure, and it is intuitively clear that the particle can 
tunnel through the barrier(s). This is manifest in the energy 
levels, because they develop a nonvanishing imaginary part 
as we vary the coupling from positive g to negative g in the 
complex plane. The smaller we choose the modulus of g, the 
bigger the humps, the longer the tunneling time of the partic- 
ular, and the smaller is the decay width of the state (i.e., the 
smaller is the modulus of the imaginary part of the resonance 
energy). Indeed, the imaginary part has a "nonperturbative" 
behaviour in the coupling and is well known to be propor- 
tional to exp(—A/\g\ b ), where A and b are positive constants. 
The same is true for the odd oscillators for g > 0; the res- 
onance energies are manifestly complex and their imaginary 
parts also involve non-analytic exponentials. 

Given the large amount of work already invested by the 
physics community into the study of anharmonic oscillators, 
it is perhaps surprising that two very basic basic questions re- 
garding the above mentioned anharmonic oscillators have not 
yet been fully addressed in the literature: (i) What are the 
higher-order corrections to the nonperturbative behaviour of 
the resonance energies, and how are the real and the imagi- 
nary part of the resonance energy described by a (possibly) 
non-analytic, generalized expansion in gl Which mathemati- 
cal structures (exponentials, logarithms, . . . ) form part of such 
an expansion? (ii) What is the general large-order behaviour 



of perturbation theory for an arbitrary energy levels of an odd 
oscillator of arbitrary degree? Both of these questions are ad- 
dressed here, and we also report on numerical calculations to 
test our analytic formulas, with special emphasis on the an- 
harmonic oscillators of degrees 6 and 7. 

Indeed, to answer the above questions, we rely in part on the 
work of Bender and Wu who, in 1971 (see Ref. (2j), solved 
question (ii) for even anharmonic oscillators, and on the con- 
cept of PT-symmetry |6j |7J HJ [9) for the formulation of a 
dispersion relation for the resonance energies of odd oscilla- 
tors. Another essential ingredient of our analysis are gener- 
alized quantization conditions which allow us to describe in- 
stanton contributions and which go beyond the ordinary Bohr- 
Sommerfeld formalism. 

Instanton actions. — We consider even anharmonic oscilla- 
tors of degree N in the convention ([TJ with energy eigen- 
values E^\g), and odd Hamiltonians in the convention ^ 
with complex resonance energies are e ( h 4 \g). Formulating 
the problem of the determination of energy levels in terms of 
a Euclidean path integral ifTOl , it becomes clear that instanton 
configurations should be analyzed, and we start with the case 
of even oscillators. Here, the instanton configuration exists 
for negative g, and we thus scale q(t) — (— sO -1 /^ -2 -* £{t). 
The Euclidean action reads 

The instanton configurations are (N even) 

q±(t) = ±(-g)-^ {1 + cosh[(7V - 2)(t - t )}}-^ . 

(4) 

Here, to is a collective coordinate. For odd anharmonic oscil- 
lators, we transform q(t) — — g _1 / ( ' 2M_4 ' ) x(£) and obtain the 
Euclidean action (M odd) 

S'[ X ] =g~^ Jdt {\x 2 + \x 2 - x M ) ■ (5) 

The instanton q(t) — q c \(t) is unique because the potential 
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FIG. 1: (color online.) Instanton world-line configuration for the 
cubic potential. The plot shows the instanton world-line Xci(t) = 
[cosh(t) + l] _1 immersed in the scaled potential U(x) = X 3 ~hx 2 - 



has lost the invariance under parity (see Fig.[T]>, 

q c i(t) = {1 + cosh[(M - 2)(t - to)]}" 1 ^ • (6) 

Inserting the solutions q^(t) and q c \(t) into Eqs. (J3J and (j5]>, 
we obtain the classical Euclidean instanton actions 

S[<£] ={-g)~^ A(N), S'[q cl ] = g-^A(M), (7a) 



2 2/( 



m-2) 



B 



m — 2 m — 2 



(7b) 



where B(x, y) — T(x) T(y) /T(x + y) is the Euler Beta func- 
tion. 

Dispersion relations. — An evaluation of the quantum fluc- 
tuations about the instanton configurations according to 
Ref. [ 10 J reveals that the imaginary part of the nth resonance 
energy En N \g), for potentials of even order N, is given by 
(in leading order, corrections are of relative order g 2 /( N ~ 2 ^) 



lmE n (N,g<0) 
x exp 
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2C{N) 



n+1/2 



n!V2yr V (- 5 ) 2/(iv - 2) 
(-g)-V(^) A(N)) , (8) 



with C(m) — 2 2 /( m_2 ). For odd potentials in the normal- 
ization given by Eq. ([2]> and for positive coupling, ignoring 
corrections of relative order g 1 /^ -2 ), 



lmei M \g>0) = 
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x exp -g 



2n\V2^ V<7 1/(M-2) 
l/(M-2) ^ (M) 



n+1/2 



(9) 



The subtracted dispersion relation [2,3, 1 1 1 for the energies 
E n (N, g) of the even anharmonic oscillators of degree N is 



1 r ds i m ^'( S +io) 

^ J-oo s(s-g) 



Using the formula ([8]> and the dispersion relation ( 10 1, one 
may calculate [2] the large-order perturbative expansion for 
the ground-state energy of the nth level of an anharmonic os- 
cillator of order N, where we use the perturbative coefficients 



in the form E^ N \g) ~ E K ^^ g K . For an arbitrary level 
of an even oscillator of arbitrary degree, we thus rederive (2) 
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which is valid up to corrections of relative order K~ x . For 
odd Hamiltonians hM{g) which involve a perturbation of the 
form s/gq M , we have only one branch cut in the energy 171 . 
The VT- symmetry for purely imaginary coupling leads to the 
following dispersion relation ll7l[T2l. 
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Based on the dispersion relation ( 12 1 and on the general re 



suit for the imaginary part of the resonance energy for an odd 
anharmonic oscillator given in Eq. ([9]), we are now in the po- 
sition to write down the large-order behaviour of the pertur- 
bative coefficients of an arbitrary level of an odd anharmonic 
oscillator of arbitrary degree. Specifically, for a resonance 
ei M \g) - Ek e( nK 9 K , we find (M odd, M > 3) 
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Subleading corrections. — In order to go beyond the 
leading-order results, further considerations are needed. 
While details of the derivation will be presented in [13 |, we 
would like to give here the essential ingredients of our for- 
malism. First of all, we scale the coordinate in the Hamil- 

0r~l _ 2[[7n/2]-m + 2 

and (|2j) as q — > g m ~ 2 q, where [X] is 

the integer part of X (m stands for a general integer). We 

then transform the Schrodinger to the Riccati equation by set- 

_ 2H?n/2]]-rra + 2 

ting f'{q)/ l f{q) — —S(q) g m ~ 2 , and we denote by 
S+(q) the component of S(q) which is even under the opera- 
tion (g, E) —> (—g, —E). For S+, we calculate the Wentzel- 
Kramers-Brioullin (WKB) expansion for S(q) for given g and 
energy E, using the algorithm described for a general poten- 
tial in Sec. 3 of Ref. 03), by expanding the solution in frac- 

2[[m/2]-m + 2 

tional powers of g while keeping the quantity g m ~ 2 E 
fixed. A recursive procedure for the construction of the WKB 
expansion of S + (z,g,E) has been outlined in Eqs. (3.40) — 
(3.42) of Ref. [ 14 j. We then integrate the WKB expansion of 
the function S + (q) around its cuts using an approach based on 
Mellin transforms as outlined in Appendix F.3 of Ref. 031 for 
general potential [the integration contour C around the cuts is 
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chosen in accordance with Eq. (3.52) of Ref. ifTH l and conjec- 
ture the following form for the result of the contour integral, 
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dzSS+{q)=A m {E,g)+\n(2^ 
lnT(^-B m (E,g))+B m (E,g) In C !J 



(14) 
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From the right-hand side, under a suitable separation of real 
and imaginary parts, one can directly read off the two func- 
tions B m (E,g) and A m (E,g), which we refer to as the 
"perturbative function" and the "instanton function," respec- 
tively. We here indicate for completeness the first few terms 
of B m (E, g) for the oscillators with M = 3, 7 and TV = 4, 6, 



B 3 (E, g)=E + g(^ + fE 2 ) + 0(g 2 ) , 
B 4 (E l9 ) = E - g (| + I E 2 ) + 0(g 2 ) , 
B 6 (E, g)=E-g(fE+lE 3 )+ 0(g 2 ) , 

R (J? n\ p I „ ( 180675 i 444381 772 

B7{E,g) -E-\-g(^ m - + -513-^ 
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(15b) 
(15c) 
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The leading term of the A functions contains the in- 
stanton action as given in Eq. ([7}: A N (E 7 g) — 
A(N) {-g)- 2/{N ~ 2) + OigW-V) for even potentials and 
A M (E,g) = A{M)g- 1 '^ M - 2 ) +0{g l ^ M - 2 ^) for odd po- 



tentials, respectively. Higher-order terms read 

A 3 (E,g) = J i- g+ g(£ + ifE 2 ) + 0(g 2 ) 
A A {E,g) = -±-g(% + ^E 2 ) + 0{g 2 ) 
A 6 (E, g) = 2 5/ 2( ^ 3) i/2 - g ( W ^ + T ^ 3 ) 

,2 f 2504899 _|_ 45769 _g-3 _|_ 17527 _g-5^ 
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The function Aj(E, g) involves the square root of the golden 
ratio 4> — (v5 + l)/2. The result for the sextic oscillator 
(N — 6) is at variance with expectation because one would 
have assumed the presence of a term of order g 2 /^ N ~ 2 ^ = 
g 1 / 2 , but it cancels, accidentally. 



In terms of the A and B functions, we can now write down 
our conjectures for the generalized quantization conditions. 
For the "stable" cases (odd oscillator, imaginary coupling and 
even oscillator, positive coupling), we simply conjecture them 
to read (we indicate two equivalent forms) 



1/ra - B m (E 7 g)]=0, B m (E 7 g)=n + 
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(17) 



where m can be even or odd. This form of the condition 
generalizes the result l/Y{\ - E) = det{H - E) = for 
the Fredholm determinant of the harmonic oscillator Hamil- 
tonian H. In the presence of instanton configurations, the 
resonance energy is slightly displaced from the energy that 
would otherwise lead to a pole of the T function according to 
Eq. ( [17] ). The displacement of E is by a nonperturbative cor- 
rection which can be evaluated exactly in leading order [see 
Eqs. ([8} and ([9}]; this means that the zero on the right-hand 
side of 1 /r[| — B m (E, g)] — has to be replaced by the non- 
perturbatively small imaginary part of the resonance energy. 
A comparison of the resulting equation to the leading terms 
of the functions A m (E, g) and B m (E, g) which emerge from 
the evaluation of the contour integral of the WKB expansion 
{14} then, in turn, naturally leads to the following conjectures 
for even potentials (g < 0) and odd potentials (g > 0), respec- 
tively, 
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In order to solve the perturbative quantization condition ( fTT} , 
we enter with an ansatz E — Eq + gE\ + g 2 E2 + . . . and 
compare coefficients in each order in g. In order to solve ([18}, 
our ansatz also has involve nonanalytic terms as implied by the 
instanton contributions, and we then expand systematically in 
powers of g and simultaneously in powers of the nonanalytic 
factor exp(— A/\g\ p ) (with A and p suitably chosen). This 
ansatz naturally leads to the triple expansions [with constant 
coefficients £j™ ,T ^- and I/ max = max(0, J — 1)], 
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The above triple expansion is involved and in need of an 
interpretation [we concentrate on the "even case" ( 19a >]. 
The term with J — recovers the basic, perturbative ex- 
pansion which has only integer powers in g. Therefore, 



(-l) K E ( n N ^. The term with J 



1 re- 



^0,0,K(N-2)/2 

covers the leading contribution in the expansion in powers of 
exp[— A(N)/(— g) 2 ^ N ~ 2 ^] to the imaginary part of the res- 
onance energy, but including perturbative corrections which 
can be expressed in terms of a fractional power series in g 
multiplying the nonanalytic exponential. The first few terms 
of this series are given below for the first excited state of 
the cubic potential. The term with J — 2 involves a loga- 
rithm of the form In ^_gy/(N-2) ) • The explicit imaginary 
parts of the logarithms cancel against the imaginary parts of 
Borel sums carried out in complex directions [16 | of pertur- 
bation series that occur in lower-order (in J) contributions. 
Indeed, the Laplace-Borel integrals have to be carried out in 
a manner consistent with the analytic continuation of the log- 
arithms ifBirrn . 

The above formulas allow, in principle, the determination 
of corrections of arbitrary order to the resonance energies of 
oscillators of arbitrary degree. We confine ourselves to only 
two results. The first result consists of the two subleading 
corrections to the first excited state of the cubic potential, 
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which have numerically large coefficients. The second result 
concerns the correction to the large-order growth of the pertur- 
bative coefficients for the ground state of the seventh-degree 
potential, 
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Note that the absence of the g 1 ^ 2 correction from the result 



( 16c 1 implies the peculiar absence of a 1/K correction to the 



leading factorial growth of perturbative coefficients for the 
sixth-degree potential (an observation which has been verified 
against high-precision numerical calculations lfT3ll ). 

Conclusions. — The general non-analytic expansions |T9]) 



for the resonance energies of even and odd oscillators are 
triple expansions in terms of non-analytic exponentials, loga- 
rithmic factors and fractional power series. They follow from 
the generalized quantization conditions (|18a|) and (18b). The 



[T8a} and 

general leading-order behaviour of perturbation theory for ar- 
bitrary levels of odd oscillators is given in Eq. (fT3|). These 



results allow us to describe the widths of arbitrary resonances 
accurately by higher-order analytic formulas. E.g., the first 



two corrections terms given in ( 20 1 are indispensable for ob- 
taining satisfactory agreement of the analytic formula for the 
first excited cubic resonance energy with numerical calcula- 
tions, even at small g ~ 0.01. 

In a wider context, the following applications of our results 
can be envisaged: in field theory (large-order estimates), the 
perturbations about the instanton configurations have usually 
been neglected in the calculation of n-point functions which 
enter the renormahzation-group equations; our results indi- 
cate that numerically large corrections may enter already on 
the level of model calculations, and it may thus be worthwhile 
to revisit the subject. Our analysis may also be helpful for the 
analytic description of resonances in quantum dot potentials 
which approximate the cubic anharmonic oscillator and are 
important for quantum computing; from our analysis, it is ob- 
vious that in a general potential, allowance should be made for 
higher-order correction terms in the description of the width 
(a functional form exp(—A/\g\ p ) [1 + B g + C g 2 + . . . ] with 
nonvanishing B and C seems to be indispensable). From a 
more fundamental point of view, we can say that in potentials 
which allow for tunneling of the particle, our analysis suggests 
that the familiar Bohr-Sommerfeld- Wilson quantization con- 
dition §pdq — 2i\nh still holds, but only if we assume a 
rather complicated analytic form for the left-hand side of this 
quantization condition [see Eq. ([14])]. 
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